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ABSTRACT 

In order to understand the orbital history of Galactic halo objects, such as globular 
clusters, authors usually assume a static potential for our Galaxy with parameters that 
appear at the present-day. According to the standard paradigm of galaxy formation, 
galaxies grow through a continuous accretion of fresh gas and a hierarchical merging 
with smaller galaxies from high redshift to the present day. This implies that the mass 
and size of disc, bulge, and halo change with time. We investigate the effect of assuming 
a live Galactic potential on the orbital history of halo objects and its consequences on 
their internal evolution. We numerically integrate backwards the equations of motion 
of different test objects located in different Galactocentric distances in both static 
and time-dependent Galactic potentials in order to see if it is possible to discriminate 
between them. We show that in a live potential, the birth of the objects, 13 Gyr 
ago, would have occurred at significantly larger Galactocentric distances, compared 
to the objects orbiting in a static potential. Based on the direct A-body calculations 
of star clusters carried out with collisional A-body code, NBODY6, we also discuss 
the consequences of the time-dependence of a Galactic potential on the early- and 
long-term evolution of star clusters in a simple way, by comparing the evolution of 
two star clusters embedded in galactic models, which represent the galaxy at present 
and 12 Gyr ago, respectively. We show that assuming a static potential over a Hubble 
time for our Galaxy as it is often done, leads to an enhancement of mass-loss, an 
overestimation of the dissolution rates of globular clusters, an underestimation of the 
final size of star clusters, and a shallower stellar mass function. 
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1 INTRODUCTION 

Tens of s atellite galax i es an d about 160 globular clus¬ 
ters IGCs: [YTar risJ l l 99(1 . l2010h , have been identified in the 
Milky Way (MW), that are distributed out to more than 
200 kpc, orbiting around the centre of our Galaxy. Ob¬ 
servations show that nearly all galaxies host these sys¬ 
tems, wi th giant ellipticals havi ng almost the largest pop¬ 
ulation (IBrodie fe Straden 120061 ). Stellar population stud- 
ies have revealed that GCs have ages up to 13 Gyr (e.g., 
IChabover fe Krausi 120021 : lHansen et al.l 120021 1 , and there¬ 
fore they represent fossil records of the earliest epoch of 
galaxy formation. As such, they are potentially powerful 
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probes of p hy sical cond itions in the high redshift Universe 
llBrodie fe Stradeill2006l ). 

All star clusters lose mass over time and this depends 
on a number of internal and external processes, as e.g., 
mass-loss due to stellar evolution, mass segregation, and 
core collapse due to two-body relaxation, and the external 
tidal field of the parent galaxy within which the cluster or- 


tidai held 01 th e parent galaxy witmn which the cluster or¬ 
bits (IVes perini &; Heggie 19971: iBaumgardt fe Makinoll2003l : 


iHeggie fe Hutl 20031 : Gieles, Heggie. fc Zhaof 201ll ). In addi¬ 

tion, the evolution of star clusters depends crucially on the 
initial conditions (e.g., the initial mass profile of the star 
cluster, the initial mass function (IMF) of the stars, and 
the initi al binary fraction) and the orbital parameters of the 
cluster j Madrid, Hurley, fe Sippell 120121 ; IWebb et al.l 120131 . 
l20l4lHaghi et al.ll2014l) 

Based on A-body simulations of GC systems, it is 
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well accepted that the GC populations we observe to¬ 
day are only the very rem n ants o f mu ch richer systems 
(e.g.. iBonaca. Geha. fc Kal ivavalill [20121: iGrillmair et al.1 


l2013i : [Brocka mp et al.ll2014l : Koposov et al.ll20l4l . The rate 


of GC erosion strongly depends on the details of the gravi¬ 
tational potential of the host galaxy as well as on the inter¬ 
nal properties of the GCs dBrockamp et al.]l2014l'l . Therefore, 
the present-day distribution of GC systems around the MW, 
and their properties may be valuable probes of the Galaxy 
potential. 

The survival or dissolution of star clusters in the galac¬ 
tic tides within which they orbit also depends crucially on 
their orbital history: star clusters with large radii spend¬ 
ing the major part of their lifetime in the innermost re¬ 
gions of our Galaxy are more susceptible to tidally in- 


ever, attempts have been made to overcome this difficulty. 
A pioneer ing study of t his s ubject was carried out by 
iRenaud. Gieles. fc Boilvl (l201lh who investigated the evo¬ 
lution of star clusters including a time-dependent poten¬ 
tial. They proposed a novel approach to extract the tidal 
information as tables of tensors from a galaxy or cosmol¬ 
ogy simulation along one orbit. Their method has been ap¬ 
plied to a large number of star c lusters in a galaxy ma¬ 
jor merger (IRenaud fc Gielesl 1201 3t l emulating the Antenna 
galaxies (NBC 4038/39). A more recent improvement of this 
method using any definition of the external potential as a 
functi on of space and time can be f ound in IRenau d fc Gielesl 
(120151) . Similar work was done bv lR.ieder et al.l ( 201314 who 


vive for a Hubble time (e.g. 

Giersz & Heggic 

1993; Hurley 

2007): Hcggie & Giersz 

2008 

Gieles. Heggie. & Zhao 201 ll: 

Brockamp et al. 2014: 

Haghi et al. 2014). Also, calculation 


of orbital motions backward in time is necessary to model 
the formation of the stellar and gaseous streams emerging 
from GCs or accreting satellite galaxies. 

Therefore, a detailed understanding of the orbital his¬ 
tory which requires a better understanding of the evolution 
of the Galactic potential since its formation, is an essential 
issue in investigating the tidal erosion of GCs (via mass- 
loss) and accreting disruption of satellite galaxies. Many 
authors usually use the static potential, i.e., assume that 
it remains unchanged during the orbital integration. But, 
observations have revealed that the size and the mass con¬ 
tent of galaxies change significantly with redshift such that 
the sizes of the galaxies at high redshifts are smaller in com- 
pari son with galaxies o f similar mass in the local universe 


tracked the tidal history of clusters in a cosmological con¬ 
text by inserting clusters into a dark matter only cold dark 
matter simulation. They also found that mergers tend to 
increase the mass-loss rates of clusters. 

In this work, we will estimate how the mass, character¬ 
istic radii, and the mass function (MF) slope of a star cluster 
would change if it is evolved in the time-dependent potential 
by comparing the evolution of two star clusters embedded 
in two galactic models, which represent our Galaxy with the 
present-day parameters and with parameter values at 12 Gyr 
ago, respectively. 

As a concrete example of motion within the Galactic 
halo, we will also use the backward motion of the Large 
Magellanic Cloud (LMC) located at about 50 kpc from the 
centre of the MW in both static and live Galactic poten¬ 
tials to see if, at least in principle, our approach is able to 
discriminate between them. 

In Section 2, we describe the characteristics of the time- 
dependent galactic potential we have used in this paper. 
We compare the orbital motions of different test particles 


(e.g 

Franx et al 

2008; Williams et al. 2010|; Moslch et al. 

around the centre of the Galaxy in Section 3. In Section 4, 

201 

; Law et al. 

20lJMosleh. Williams,& Franx 2013!). 

we present results of A-body calculations for the dynamical 


There are 

many proposed scenarios to explain the 

evolution of star clusters moving through an external galaxy 


physical processes of galaxy assembly that well reproduce 
the observable properties like, e.g., the stellar mass and 
size of galaxies at different redshifts. Among them are the 


galaxy minor or major mergers (e.g., Khochfar&^jilk[ 2009; 
iKhochfar fc Silkll20()9l : iNaab. Johansson, fc Qstrikeill2009l) . 
and the accretion of fresh gas i n outer regions activati ng new 


star formation (e.g., Elmegreen, Bo urnaud. fc Elmegreenl 
120081 : iDekel, Sari, fc Ceverinol 120091 ). Indeed. in the stan¬ 
dard picture of galaxy formation, galaxies are em¬ 
bedded in massive virialized haloes of dark matter 
(ISpringel. Frenk, fc Whitell2006l ). These dark matter haloes 
accumulate over time hierarchically, continuously growing 
via accretion of dark matter and merging with other haloes 
from high redshift to the present day. The fraction of GCs 
that have survived to the present day has evolved in a time- 
dependent potential of the host galaxy. 

Understanding the influence of the time-dependence of 
the Galactic potential on the orbital history of the halo ob¬ 
jects (e.g., GCs) rotating around the MW at different Galac- 
tocentric distances, and its consequence on their early- and 
long-term evolution over a Hubble time is the main motiva¬ 
tion for this paper. 

Indeed, because of the wide range of time-scales and 
size scales, from the two-body strong encounters of stars 
to the galactic scales, the iV-body simulations of star 
clusters in a live galactic potential are challenged. How- 


with different background potential parameters. The simu¬ 
lations were carried out with the collisional N-body code 
nbodyG on desktop workstations with Nvidia 690 Graphics 
Processing LTnits at the Institute for Advanced Studies in 
Basic Sciences (IASBS). Finally, in Section 5, we summa¬ 
rize our results. 


2 THE GALACTIC POTENTIAL 
2.1 Static Galactic potential 

In the first part, for the Galaxy potential we, like many 
other authors, assumed that the galaxy potential is static 
and consists of three idealized components, 


$tot = + $6 + 


( 1 ) 


including _ a Miam oto-Nagai disc potential 

llMivamoto fc Nagnilll975l ) given by 


<S> d (x,y,z) = - 


GM d 


\Jx 2 + y 2 + (a + %/z 2 + 


a central bulge, the Hernquist (1990) model, given by 

GMb 


$ b{x,y,z ) = - 


r + r c 


( 2 ) 


( 3 ) 
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Figure 1. Left-hand panel: logarithmic plot of the size scale of the MW’s mass components as a function of cosmological time. The virial 
radius of the dark matter halo r v i r (dashed green), scale radius of the disc a (dotted blue), the scale-height of the disc b (dash-dotted 
dark blue), and the scale radius of the bulge r c (solid red) are shown in the figure. Right-hand Panel: the evolution of the mass scale of 
the MW’s mass components as a function of cosmological time. The virial mass of the dark matter halo M v i r (dashed green), the scale 
mass of the disc (dotted blue), and of the bulge M & (solid red) are shown in the figure. All components show fast growth in the past 
time. 


Disc 

Bulge 

Halo 

M d = 7.5 x 10 10 

M b = 2.5 x 10 10 

M v i r = 9 X 10 11 

a = 5.4; b = 0.3 

r c = 0.5 

r v i r = 250 



c = 13.1 


Table 1. The present-day parameters of the mass components 
of the MW-like potential used in our calculations. Masses and 
distances are in Mq and in kpc, respectively. 


Disc 

Bulge 

Halo 

M d = 6.6 x 10 9 

M b = 2.2 x 10 9 

M v i r = 7.9 x 10 10 

a = 0.65; b = 0.04 

r c = 0.06 

r v i T = 29.9 



c = 2.86 


Table 2. The parameters of the mass components of the MW-like 
potential at T = —12 Gyr. Masses and distances are in Mq and 
in kpc, respectively. Eqs. 5-9 are used to calculate these values 
at z = 3.5 (i.e., T = —12 Gyr). 


and a NFW dark matter halo (| Navarro, Frenk, &; Whit el 
EssB) with a potential of the form 


4’iV.FW 


GMvir 

r[log(l + c)- jqy 



(4) 


Here, r = x 1 + y 2 + z 2 is the distance from the galactic 
centre at any given time; Mb, Md and M v i r are the char¬ 
acteristic mass of the bulge, disc, and halo, respectively; r c 
and r v i r are the characteristic radius of the bulge and halo, 
respectively; and a is the scale radius and b the scale height 
that adjust the shape of the disc. The present-day numerical 
values of these parameters at redshift z = 0 are given in the 
Table Q] 


2.2 Time-dependent Galactic potential 


To model the evolution of the Galactic potential, we as¬ 
sume a semicosmological time-dependent gravitational po¬ 
tential in which the characteristic parameters vary in time. 
The evolution of the mass and the virial concentration of 
the galaxy’s halo as a function of redshift are given b y 
llWechsler et al.ll2002l : IZhao et aPl2003l : iGomez et al.ll201oh 

M vir (z) = M vir (0) exp(—2a c 2), (5) 

where the formation epoch is set to be a c = 0.34, and 
c(0) 


c(z) = 


1 + z' 


( 6 ) 


For the disc and b ulge we follow the recipe given by 
iBullock fe Johnstonl (l2005l 'l. as for masses we have 


Md,b{z) = M v i r (z) 


Md,b(0) 


(7) 


(8) 


Mvir [0) ’ 

and for scale-lengths the evolution can be expressed as 
{a, b, r c }(z) = r v i r (z) ^ ^ . 

i viryyJ) 

Here r V i r is the virial radius of the dark matter halo, varying 
as 

r. (») 


47tA V ir{z)pc{z) J 
where A v i r (z) denotes the virial overdensity , 
A v i r (z) = 18t r 2 + 82[fiQ) - 1] - 39[fiQ) - l] 2 
with fiQ) the mass density of the universe, 
flm.o (1 + z ) 3 


( 10 ) 


( 11 ) 


^ ^ O m ,o(l + z) 3 + 11a, o ’ 

and p c {z) is the critical density of the universe at a given 
redshift, 


pc(z) = 


with 


3 H 2 {z) 
87rG 


( 12 ) 
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Figure 2. Comparison of the circular velocity curve as a function 
of Galactocentric distance at 2 : = 0 (black solid line) with rotation 
curve at z = 3.5 (grey dashed line). The asymptotic values of 
circular velocity in two frame work differ by about 75 kms -1 
after about 12 Gyr of evolution. 


H{z) = Ho v^A.o + tWtl + z) 3 . (13) 

We also adopted a flat cosmology defined by 
= 0.3 and 17 a, o = 0.7 with a Hubble constant 
of H(z = 0) = H 0 = 70kms” 1 Mpc” 1 . 


In cosmology one can label the time t since the big bang 
in terms of the redshift of light emitted at t. By integrating 
the Friedmann equation, the behaviour of the cosmological 
redshift in terms of time for a flat universe can be found as 


- 1/3 


/ Sim, o sinh 2 (| Hot\J 17a,o) 
\ SI A,0 


(14) 
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In Fig. [2 we display how the characteristic parameters 
of the MW (i.e., all Galactic mass components and scale- 
lengths) vary as a function of cosmological time. The pa¬ 
rameters of the mass components of the MW-like potential 
at T = —12 Gyr (z = 3.5) are given in Table 2. 

The Galactic rotation curves deduced from the three- 
component mass model, described in Eqs. 2-13, at z = 0 and 
z = 3.5 are plotted in Fig. [5] As can be seen, the present- 
day circular velocity curve of this model takes a value of 
about 221 kms” 1 at 8.5 kpc from the galactic centre, while 
it is about 117 kms -1 at z = 3.5. The smaller value of 
asymptotic rotational velocity at Rg = 100 kpc occurs in 
the time-dependent gravitational potential, while the static 
potential yields the larger one, differing by about 75 kms” 1 
after 12 Gyr of evolution. 


3 MOTIONS BACKWARD IN TIME 

In this section, we describe the results from the numerical in¬ 
tegration of the equation of motion to find the trajectory of a 
test object located in different Galactocentric distances mov¬ 
ing under the static Galactic potential and compare them 
with the same calculations in a time-dependent Galactic po¬ 
tential by using Eqs. 5-13 and the present-day parameters 
given in Table 1. 

Using the gravitational potential components described 


Figure 3. Upper panel: planar shape of the orbits of a test par¬ 
ticle w.r.t. the MW, starting from Re = 4 kpc (the first line of 
Table[3]l . The red dashed line shows the trajectory of test object 
within the time-dependent potential while the blue line is for the 
static potential. Lower panel: the time evolution (for the past 13 
Gyr) of the radial distance to the Galactic centre of a test par¬ 
ticle, starting from Rq = 4 kpc in the time-dependent potential 
(red dashed line) and the static potential (blue line). 


above, we will numerically integrate the three scalar differen¬ 
tial equations, written in Cartesian coordinates, correspond¬ 
ing to the vector differential equation, 

? = —V€>tot- (15) 

Here we use the fourth-order Runge-Kutta method to cal¬ 
culate the equation of motion of test objects and then to 
extract the trajectory of object by backtracking orbit from 
its current position and velocity. Using the initial conditions 
of Tabic [3] the equations of motion are then numerically in¬ 
tegrated backward in time for 13 Gyr in both live and static 
gravitational potentials. In this manner, we therefore obtain 
sets of initial positions and velocities required for forward 
integration in time and can be applied in e.g., A-body sim¬ 
ulations of realistic GCs or satellite galaxies. 

We calculate the orbit for three test objects: an inner 
test object currently located at Rg = 4 kpc, a solar-distance 
object located at Rg =8.5 kpc, and an outer object with 
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Figure 4. The same as Fig. [3] but for a test object that located 
initially at Rq = 8.5 kpc as shown in the second line of Table [3] 
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the present-day Galactocentric distance of Rg = 100 kpc. 
For all objects, we calculate the trajectory within two dif¬ 
ferent Galactic models: the static Galactic potential and 
time-dependent Galactic potential. The orbital parameters 
of these test objects are summarized in Table [3] 

The initial velocities are extracted from the present-day 
rotation curve. That is, the orbits are circular in the static 
potential, while in the live potential the Galactocentric dis¬ 
tances increase looking backwards in time. Note that, the 
test objects evolve on circular orbits at different galactocen¬ 
tric distances in the disc plane (i.e., the inclination angle 
of the orbits w.r.t the galactic disc is zero). In Figs. [3]-[S] 
we plot the orbital sections in the Galactocentric coordi¬ 
nate planes of the numerically integrated trajectories of the 
test objects from now to 13 Gyr ago. The time evolution of 
distances from the centre of the MW for all test objects are 
shown in the lower panel of Figs. [3]-[5] As can be seen in Fig. 
[S] the Galactocentric distances of these objects at t = —13 
Gyr are larger than initial values at t = 0 by a factor of 
about 2 (for test object initially located at Rg = 100 kpc), 
and 2.5 (for test objects initially located at 4 and 8.5 kpc). 
It should be noted such a large difference in determining the 
birth place of realistic objects (like e.g., GCs in the MW) 



-200 -100 0 100 200 
X [Kpc] 



Figure 5. The same as Fig. [3] but for a test object initially 
located at Rg =100 kpc as shown in the third line of Table[3] 



t [Gyr] 

Figure 6. The evolution of the relative difference of the Galacto¬ 
centric distance in the static and live Galactic potential in time 
that are plotted for individual test objects moving in different 
initial distances from the centre of Galaxy as given in Table [3] 
The differences are small at the beginning of orbital backward 
motion and starts to increase at about t = —6 Gyr. The value of 
Galactocentric distance get doubled within a Hubble time for a 
test object initially located at Rq = 100 kpc, while it increases 
by a factor of about 2.5 for other two test objects at Rq = 4 and 
8.5 kpc. 


may pose problems concerning the dynamical evolution. I 
will be back to this important issue in more details in Sec. 

m 


3.1 Motion of the LMC backward in time 

As a concrete example of motion within the Galactic 
halo, let us consider the orbital motion of the LMC 






























6 Haghi et al. 


Table 3. Initial coordinates, in kpc, and initial velocities, in 
kms —1 , in Galactocentric rest frame adopted for different test 
objects used in this paper. We have used the present-day rota¬ 
tion curve of our Galaxy for initial circular velocities. 



r(x,y,z) 

3D v(x, y , z) 


model 1 (inner-part object) 

(4,0,0) 

(0,215,0) 

77 

Oh 

model2 (solar-distance) 

(8.5,0,0) 

(0,221,0) 

W 

model3 (outer-part object) 

(100,0,0) 

(0,195,0) 

N 


Table 4. The coordinate, in kpc, and velocity component, in km 
s _1 , of LMC used in this work in a Galactocentric rest frame with 
the z-axis pointing towards the North Galactic Pole (NGP), the x- 
axis pointing in the direction from the Sun to the Galactic centre, 
and the positive r/-axis is directe d towards the Sun’s Galactic 
rotation l lMastropietro et al.ll2005lh 


200 


100 


-100 


-200 
-200 


V/v' ff 

/ 

✓ A f 

/ /\ f 

■ Live-Potential 

Static-Potential 

\ 

\ 

'NiA \ l 

I A' / / 

/ / \ 31 

/ / \i\ 

/ / NLX 

i 1 

* l A \ -— 

i \ 'a' 

1 \ \ 

\ v'v-' 

/ 1 / 

iA/ 

— \ i 

// ! 

7 i 
// 

/ / 


-100 


100 


200 


r(x,y,z ) 3D »(i, j, z) 

(0,-43.9,-25.04) (-4.3,-182.45,169.8) 


currently located at about 50 kpc from the centre of the 
MW. Many authors have shown that the position of the 
Magell anic Stream (M S) Q follows the orbits of M Cs ( see 


e.g. 

iMoore & Davis 19941: Connors. Kawata. & Gibson] 

2006 

: lHaghi, Rahvar. & Zonoozi 2006; Besla et al.! 

2007 

: Haghi. Zonoozi. & Rahvar 20091: Haghi & Rahvail 

201(j 

MS 

). Moreover, the shape and the kinematics of the 
is strongly influenced by the overall properties 


of the underlying potential llMurai & Fuiimoto 

198C 

Lin & Lynden-Bell 1982i; Heller & Rohlfs 19941: Sofue 

1994; 

Gardiner & Noguchi 199(1). Therefore, it is most striking to 


compare the trajectory of the LMC with orbits predicted in 
both static and live Galactic potentials. Table 2] summarizes 
the present-day Galac tocentric Cartesian coor dinates and 
velocities of the LMC llMastropietro et alj|2005l) . 

Applying the same method for the LMC to numerically 
integrate the equation of motion, we extract its trajectory 
by backtracking orbit from its current position and velocity 
for 13 Gyr. Fig. |T] depicts the planar shape of the orbits of 
the LMC in Y — Z plane (upper panel). We also show the 
evolution of the distance to the MW of the LMC in both 
static and live Galactic potentials in the lower panel of Fig. 
□ 

Another issue which should, in principle, take into ac¬ 
count is the number of disc (i.e., the MW’s disc) passage in 
both Galactic potentials. In fact, crossing the disc would im¬ 
ply a strong orbital perturbation of the LMC and perhaps 
a gas shock which can lead to the formation of stars. We 
found that the number of disc passage with Galactocentric 
distances smaller than 100 kpc, in the static potential is five, 
while it is three in the live potential. 


1 A narrow band of neutral hydrogen clouds lies along a great 
circle from (Z = 91°; b = —40°) to (Z = 299°; b = —70°), started 
from the MCs and oriented towards the South Galactic Pole. 


Y [Kpc] 



Figure 7. Top: section in the Z — Y plane of the integrated tra¬ 
jectories of the LMC for the static potential (blue line) and the 
live potential (red dashed line). The time span of the integra¬ 
tion is —13 Gyr ^ t ^ 0. The initial conditions fo r position and 
velocity are taken from (Mastropictro ct al. 200j]). Bottom: the 
evolution of the Galactocentric distance of LMC as a function 
of time. The red dashed line: time-dependent Galactic potential. 
Blue line: static potential. The galactocentric distance of LMC, 
13 Gyr ago, would have been at Rq = 60 kpc for Galactic model 
with static potential, while in the time-dependent Galactic poten¬ 
tial it would have been at larger Galactocentric distance of about 
Rq = 230 kpc. 


4 DYNAMICAL EVOLUTION OF STAR 

CLUSTERS 

Our Galaxy hosts around 160 GCs. It is well understood 
that the gravitational potential of host galaxy has a direct 
influence on the survival and the evolution of GCs: they lose 
stars through tidal stri pping and disc shocking. As already 
have shown bv lPraagman, Hurley, fc Power] ]20ld ). varying 
the mass and concentration of the halo affects the rate at 
which the star cluster loses mass. Using A-body models of 
low-number star clusters, they found that increasing the halo 
mass and concentration drives enhanced mass-loss rates and, 
in principle, implies shorter dissolution time-scales. 

Several studies have been addressed the evolution 
of clusters in time-dependent galactic potential by ar¬ 
bitrary switching tidal effects to mimic the accretion 

of a dwarf satellite on to a massive ho st _ galaxy 

llMiholics. Webb, fe Sills! I20l4 iBianchini et al.l 120151) . or 
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Figure 8. The evolution of total mass (left-hand panel) and half-mass radius (right-hand panel) with time for two simulated star 
clusters orbiting within different Galactic potentials with parameters listed in Tables [l] and [2] Cluster moving in a circular orbit with a 
Galactocentric radius of Rg = 15.5 kpc and affected by the gravitational field of our Galaxy with parameters at z = 3.5 (corresponding 
to T = —12 Gyr) is shown by a black dashed line, and with Galactocentric radius of Rq = 8.5 kpc and affected by the gravitational 
field of our galaxy with parameters at z = 0 (corresponding to the present day) is shown by a red solid line. For the cluster close to the 
galactic centre, expansion is limited by the strong tidal field and dissolves before a Hubble time. 


by rapidly varying evolutions (e.g., galaxy interac¬ 
tions and mergers) combi ning the galaxy simulati o ns to 
star cluster simulations (|Renaud, Gieles, fc Boilvl l201ll : 
iRenaud fc Gieled 120131 . 1 20151 ), or by evo lvin g sta r clusters 
in a cosmological environment (IRieder et all 12013 1. 

This issue has been simplified by several authors. For 
example. iMadrid. Hurley, fc Martial (J201J) have studied the 
impact of the host disc mass and geometry on the survival of 
star clusters by means of iV-body simulations. They showed 
that a more massive disc enhances the mass-loss rate of an 
orbiting star cluster owing to a stronger tidal field such that 
doubling the mass of the disc halves the dissolution time of a 
star cluster located at Rg = 6 kpc from the centre of Galaxy. 
They placed several of these simulations together, each time 
increasing the galaxy’s mass, to represent a realistic mass 
growth history of the MW driven by mergers of satellite 
galaxies. 

As we have shown by backward tracking of test ob¬ 
jects, in the frame work of time-varying Galactic potential, 
the Galactocentric distance of a test object (like e.g., a star 
cluster in the MW) at t — —12 Gyr is quite a bit larger than 
initial values at t = 0, by a factor of about 2. In other words, 
in order to be a star cluster in its current position and veloc¬ 
ity, it should be located at a larger Galactocentric distance 
in the past, i.e., 12 Gyr ago, when the Galaxy was also much 
lighter than its present-day mass. This implies that the clus¬ 
ter presumably was tidally underfilling 0 at the beginning of 
its evolution. Therefore, the slower mass-loss rate of clusters 
initially lying inside their tidal radii, takes a longer time to 
lose a given amou nt of mass in comparison to tidally filling 
clusters (see e.g ., iBaumgardt, De Marchi, fc Kroupal 120081 : 
iMarks. Kroupa. fe Baumgardt l2008h . Indeed, in tidally lim¬ 
ited clusters, the early evolution of massive stars leads to a 
rapid expansion, and hence a larger flow of mass over the 


2 with rh/rt values smaller than 0.05, where ry, and rt are the 
cluster half-mass and tidal radii, respectively. 


tidal boundary. It may therefore help to dissolve them more 
rapidly. Tidally underfilling clusters, however, can survive 
this early expansion. 

We therefore expect that different scenarios for Galac¬ 
tic potential (i.e., the time-varying versus the static invari¬ 
ant Galactic potentials), in principle, can lead to different 
evolution and survival of star clusters and consequently dif¬ 
ferent depletion rates of satellite star clusters. Here, in this 
section we assess this difference by direct A r -body simula¬ 
tions of GCs in a realistic MW-like potential using the code 
nbody6. 

The current version of nbody6 does not allow for a 
treatment of tides with an explicit time-dependent back¬ 
ground potential. Within the current framework of nbody6 
and in order to estimate the fate of a star cluster which its 
host galaxy grows with time we have calculated two inde¬ 
pendent models of star clusters with different masses and 
sizes for host galaxy. 

• First we assume a star cluster in a circular orbit with 
a Galactocentric distance of Rg = 8.5 kpc that evolves a 
’’heavy” galaxy with a mass and geometrical parameters of 
the present-day MW-like potential listed in Table 1. 

• We then simulate the same star cluster in a circular or¬ 
bit with Galactocentric distance of Rg = 15.5 kpc within a 
"light” host galaxy with parameters of the MW-like poten¬ 
tial at T = —12Gyr given in Tabled The reason for this 
choice of Galactocentric distance is that a cluster which is 
currently located at 8.5 kpc from the centre of the Galaxy 
had previously started to evolve (12 Gyr ago) at a Galacto¬ 
centric distance of Rg = 15.5 kpc within the time-varying 
Galactic potential (See Fig. [4]). 

Our simulations include mass-loss driven by stellar evo¬ 
lution by using the SSE/ BSE rout i nes an d anal ytica l fit¬ 
ting functions developed by iHurlev, Pols fe Toutl (l2000l l and 
IHurlev, Tout fe Pols! (120021 1. the two-body relaxation, and a 
realistic treatment of the external tidal field. Both star clus¬ 
ters are evolving with initial particle number of IV = 50000 
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Figure 9. The evolution of the core (lower curves) and tidal ra¬ 
dius (upper curves) is plotted as a function of time in logarithmic 
scale. Models are the same as Fig. [8] The tidal radius of cluster 
with orbit at Rq = 8.5 kpc from the galactic centre decreases 
continuously due to ongoing mass-loss. Core collapse is reached 
at T ~ 9 Gyr for this cluster. But, the star cluster that evolves in 
a "light” galaxy on a circular orbit with radius of Rq = 15.5 kpc 
does not exhibit core collapse before a Hubble time. 


(corresponding to M ^ 30000 Mq ) that were distributed as 
a Plummer density profile (Plummer 1911) with the same 
initial half-mass radius of r^.p = 3 pc. The models started 
with a Kroupa stellar IMF (jKroupal 1200 ll : iKroupa et al .1 
l2013h . which consists of two power laws with slope Qi = 1.3 
for stars with masses between 0.08 and 0.5 Mq and slope 
02 = 2.3 for more massive stars. The range of stellar masses 
was chosen to be from 0.08 to 100 Mq. The simulated clus¬ 
ters evolve on circular orbits at different galactocentric dis¬ 
tances in the disc plane (i.e., the inclination angle of the 
orbits with respect to the galactic disc is zero). Our main 
focus is to study how quantities which can be checked ob- 
servationally, as e.g., the slope of the MF or the size scale 
of the cluster change with time. 

4.1 The evolution of mass and characteristic radii 

The evolution of the total mass of the simulated clusters 
with time for two different Galactic parameters as given in 
Tables |T| and [2] are plotted in Fig. [8] In fact, the long-term 
mass-loss for these clusters can be regarded as a runaway 
overflow over the tidal boundary. It can be seen, a star clus¬ 
ter that evolves in a "light” galaxy on a circular orbit with 
radius of Rg = 15.5 kpc has remaining mass of 10 4 Mq, 
i.e., 30% of its initial mass after a Hubble time of evolu¬ 
tion, while a star cluster evolving in a ’’heavy” galaxy with 
Rg = 8.5 kpc does not last for a Hubble time and dissolved 
after about 11 Gyr. 

Fig. [8] (right-hand panel) depicts the time-evolution of 
the 3D half-mass radius of simulated star clusters at different 
galactocentric distances. Clusters start with an initial half¬ 
mass radius of 3 pc and undergo an expansion triggered by 
stellar evolution within the first Gyr. The half-mass radius of 
a simulated star cluster evolving within a ’’heavy” Galaxy 
(with the present-day values of parameters of MW listed 
in Table [TJ at 8.5 kpc from the galactic centre, reaches a 



Figure 10. Shown are the evolution of Lagrange radii and core 
radius of a cluster that evolves in a ’’heavy” galaxy with orbit 
at Rg = 8.5 kpc from the galactic centre as a function of time. 
LR2, LR5, LR10, and LR50 are, respectively, the 2, 5. 10 and 50 
per cent Lagrange radii. Both the half-mass and core radii show 
an initial increase corresponding to stellar evolution and mass-loss 
from massive stars. Cluster reaches the end of initial core-collapse 
phase at T ~ 9 Gyr. The sharp change in the behaviour of the core 
at T ~ 9 Gyr, when the size of the core suddenly increased might 
be linked to an interaction within the core involving a dynamically 
formed binary. 


maximum value, which appears to be clearly linked to its 
galactocentric distance, before it decreases again until the 
cluster dissolves after 11 Gyr of evolution. As shown in Fig. 
[8] the evolution of a star cluster orbiting at Rg =15.5 kpc 
under the tidal field of a ” light” Galaxy with the parameters 
of the MW-like potential at T = —12 Gyr (Tabic [lj appears 
to have not reached its tidal limit yet and keeps expanding 
till the end of the simulation. This is easy to understand: 
a smaller orbital radius leads to a faster disruption and a 
smaller half-mass radius after 13 Gyr evolution owing to 
the enhanced mass-loss driven by the galactic tide, and the 
stronger cut-off it inflicts on the clusters. 

Another useful diagnostic of the difference between 
models at different galactocentric distances is provided by 
the core 0 and tidal radii. Their time evolutions arc illus¬ 
trated in Fig. [U] Accentuated mass-loss at Rg = 8.5 kpc as 
seen in Fig. [8] precipitates the onset of core-collapse. The 
time of core-collapse is usually determined by the moment 
of core bounce, which is seen in the tim e ev olu tion ofjthe 
core radius or density radius (iFuiii fc Porteeies Zwartll2014l f. 
However for some models it is difficult to distinguish the 
core-collapse, because there does not seem to be a peak in 
the density evolution or a depression in the core radius. The 
core radius shows a minimum at ~ 9 Gyr for a star cluster 
evolves at Rg = 8.5 kpc from the centre of ’’heavy” galaxy 
model which we identify as the moment of the initial core- 

3 In our simulations, the core radius is a density weighted aver¬ 
age distance of ea ch star to the point of highest stellar density 
within the cluster iCasertano fe Hutlll985l : lAarse thl l2003h , while 
in observational studies the core radius is generally defined as the 
radius where the surface brightness falls to half its central value 
(King 1962). 
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collapse phase ends (Fig. [9j). According to Fig. 1101 the core 
radius expands in the beginning (within the first 2 Gyr) due 
to the weak tidal field, but eventually at T ~ 9 Gyr it goes 
into a small core collapse. But since it is also close to be¬ 
ing disrupted, the core collapse does not look like the usual 
deep contraction. The cluster expands in a jump, probably 
due to the ejection of the central binary or even a binary¬ 
binary ejection. This leaves the core much less bound than 
before the ejection, and hence the cluster expands and dis¬ 
solves in an instant. The evolution of Lagrange radii also 
supports this conclusion as they go up at T ~ 9 Gyr, and 
one can conclude this is core collapse. Short-term effects on 
the evolution of core radius are often related to the presence 
of high-energetic binaries, probably comprised of two black 
holes in the core. 

The cluster in orbit of 15.5 kpc (” light” Galaxy) does 
not reach the end of the core-collapse phase in 13 Gyr. 
Therefore, we conclude that the core radius is affected 
by c hang es in tidal forces on the cluster, in agreement 
with [Madrid, Hurley, fe Sippell d2012h who have found that 
the galactocentric distance of a star cluster has an im¬ 
pact on its core radius and the onset of core collapse. 
This is in contradictio n with the conclusion reached by 
iMiholics. Webb, fe Sills! d2014l ! who showed that the core 
radius of a cluster will depend on the initial structural con¬ 
ditions of the cluster and will not be affected by its tidal 
history. It can be seen that the tidal radius of the cluster 
in the heavier (present day) galaxy is smaller (tidally lim¬ 
ited) than the cluster in the lighter (12 Gyr ago) galaxy. We 
therefore conclude that clusters survive longer in an evolving 
galaxy than in a galaxy which is kept static. 

As expected, galaxy with masses between these two ex¬ 
tremes (i.e., a time-dependent galaxy mass model) define in¬ 
termediate regime of mass-loss. One can therefore conclude 
that after a Hubble time of evolution, a star cluster has less 
mass in a static potential than a simulated cluster evolving 
in a live potential. This is because, in the static potential 
the galaxy mass on average (over time) is larger than the 
live potential. As already have shown by several authors, 
increasing the mass of h ost galaxy accelerate the de struc¬ 
tion time of star clusters dPraagman, Hurley, fe Powei|[201(] : 
iMadrid. Hurley, fe Martidl2014| j. 

However, the evolution of a star cluster strongly de¬ 
pends on its filling factor ( rn/rt )• The present-day under¬ 
filling clusters (i.e., rh/rt < 0.05) would remain underfilling 
in an evolving-potential over the whole 12 Gyr of evolution. 
This is because of the larger Rg and lighter Mg compared 
to a static potential (see Fig. 0. So one would expect the 
cluster’s final size to be nearly similar to that of the static 
case. But, a present-day tidally limited cluster would be un¬ 
derfilling for a while in an evolving galactic potential. 

4.2 The evolution of the MF 

The stellar MF is one of the most important observable 
parameters that changes through the dynamical evolution 
of star clusters. It is evident that the two-body relaxation 
driven evaporation through the tidal boundary of a star clus¬ 
ter gives rise to a significant correlation between the MF- 
slope and the strength of tidal field of host galaxy. Because 
of the dynamical mass segregation occurs due to two-body 
relaxation, the evaporation rate is larger for low-mass stars 



Time [Gyr] 



Figure 11. The evolution of the global stellar MF-slope for low- 
mass (m C 0.5 Mq, top panel) and high-mass stars (m ^ 0.5M@, 
bottom panel) within the half-mass radius is plotted as a function 
of time. Models are the same as Fig. [8] A flatter MF is the direct 
result of increased tidal stripping of outer region stars that are 
preferentially low in mass due to dynamical mass segregation. 


than it is for high-mass ones (Giersz & Heggie 1997, Baum- 
gardt & Makino 2003). Thus, the preferential escape of low- 
mass stars leads to the flattening of the MF as the dynamical 
evolution of a star cluster proceeds. It is shown that the evo¬ 
lution of MF-slope is faster for clusters at smaller galacto¬ 
centric distances (Vesperini & Heggie 1997), i.e., experience 
a stronger tidal field (Webb et al. 2013). 

The canonical IMF as observed in young star clusters in 
the MW is often expressed as a two-part power-law function 
(<fm ^ m ~ a ) with near Salpeter-like slope above O.5M0 
(i.e., Q 2 = 2.3: ISal pete till 955)l i. and a shallower slope of an = 
1.3 for stars in the mass range 0.08 — 0.5 M© [|Kroupall200ll : 
iKrouna et al .1 [jobl ). 

Fig-[H]depicts the evolution of the MF-slope at the low- 
mass end as a function of time. The slope was determined 
from a fit to the distribution of stars with masses m ^ 0.5. In 
both cases the MF flattens as the cluster loses stars. There¬ 
fore, as the galactocentric distance of a GC decreases, the 
strength of the stellar mass-loss driven by two-body relax¬ 
ation increases, and hence the amount of the flattening of the 
MF enhances. Fig. II II confirms that the slope of the MF in 
a cluster orbiting at Rg = 8.5 kpc within a ’’heavy” Galaxy 
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changes significantly as compared to a cluster evolving in a 
circular orbit with a Galactocentric radius of Rg = 15.5 kpc 
within a ” light” galaxy model. 


Recent observational work on a number of MW GCs 
have shown that the global MF-slope in the low-mass range 
is significantly shallower than a canonical MF-slope of about 


2.3 (Kroupa 2001, see e.g. 

Dc Marchi et al. 2007: Jordi et al. 

200£ 

: Paust et al. 2010l: 

Frank et al. 20121: Hamren et al. 

2013 

). However, the preferential loss of low-mass stars due 


to two-body relaxation would b e a natural explanation fo r 
the observed MF depletion llBaumgardt fe Makinol [20031! . 
for diffuse outer halo clusters such as Pal 4 and Pal 14 
(i.e., a low mass together with a large half-mass radius), 
the present-day two-body relaxation time is of the order of 
a Hubble time. Therefore, relaxation should be inefficient in 
these clusters and the observati ons should be an indicatio n 
for primordial mass segregation dZonoozi et al HB. 
Our findings in this section show that considering a live po¬ 
tential for our Galaxy makes it even more difficult to explain 
the observed MF flattening. This is because, the changes in 
the MF due to the tidal stripping is less in a live potential 
compared to the static potential model. 


5 CONCLUSIONS 

Many authors usually use the static potential for our Galaxy, 
as the common assumption is that it remains unchanged 
during the orbital integration. In this paper, we have inves¬ 
tigated the influence of the time-dependence of the Galactic 
potential on the orbital history of the halo objects and its 
consequences on their internal evolution. First, we numer¬ 
ically integrated backwards the orbits of different test ob¬ 
jects over a Hubble time, located in different Galactocentric 
distances within both static and live (cosmologically moti¬ 
vated) Galactic potentials to assess the possible differences. 

It turns out that, the static and live potential do yield 
different trajectories for our test objects orbiting in differ¬ 
ent Galactocentric distances. We have shown that the spa¬ 
tial extinction of the orbit’s section in the coordinate planes 
is larger for live potential w.r.t static potential, such that 
in a live potential, the birth of the objects, 13 Gyr ago, 
would have occurred at significantly larger Galactocentric 
distances, compared to the objects orbiting in a static po¬ 
tential. 

As a concrete example of motion within the Galactic 
halo, we also used the backward motion of the LMC in both 
static and live Galactic potentials. In addition to the differ¬ 
ent trajectories of LMC we uncovered here, we found that 
the orbital period of the LMC around the MW is about 2.7 
Gyr in the static potential model, while it is 4 Gyr in a 
time-dependent Galactic potential model. This is important 
because it is believed that the kinematics and morphology 
of the MS follows the orbits of the LMC. We furthermore 
found that, in the static potential, the number of disc pas¬ 
sage with Galactocentric distances smaller than 100 kpc is 
five, while it is three in the live potential. 

We finally investigated the impact of the assuming a 
live potential on the dynamical evolution of star clusters 
by means of the collisional TV-body code, nbody6. Since 
the current version of Nbody6 does not allow for a treat¬ 
ment of live potential of host galaxy with an explicit time- 


dependence, we calculated two models of star cluster with 
different masses and sizes for host galaxy (which represents 
the galaxy at present and 12 Gyr ago), to roughly estimate 
how a cluster’s half-mass radius, the total mass, and the 
MF-slope develops over the time. 

We followed the evolution of clusters at different Galac¬ 
tocentric distances with different Galaxy mass models, and 
found that the weaker mass-loss of clusters evolving in a 
weaker tidal field (i.e., at a larger Galactocentric distance 
and within a light-mass Galaxy) leads to a significantly 
larger final size. 

Our computations demonstrate that for two star clus¬ 
ters moving in circular orbits with different Galactocen¬ 
tric distances and within different Galactic models, one 
with Rg = 8.5 kpc and within a heavy-mass Galaxy (with 
the present-day parameters of our Galaxy) and one with 
Rg = 15.5 kpc evolving in a light-mass Galaxy (with pa¬ 
rameter values at 12 Gyr ago), the star cluster with smaller 
Galactocentric distance will have: 

• a stronger tidal truncation and a smaller size 

• an enhanced mass-loss rate and a shorter dissolution 
time 

• a flatter MF 


Since a galaxy with parameters between these two ex¬ 
tremes defines intermediate regime of mass-loss, one can 
therefore conclude that over a Hubble time of evolution in 
a semicosmological time-dependent Galaxy model, the star 
cluster has more mass than a simulated cluster evolving in a 
static invariant potential; this is because in the static poten¬ 
tial the galaxy mass on average (over time) is larger than the 
live potential. This implies that assuming a static potential 
for our Galaxy (as it is often done) leads to an enhance¬ 
ment of mass-loss rate, an overestimation of the dissolution 
rates of GCs, and an underestimation of the final size of star 
clusters. 

Consequently, after a Hubble time of evolution in the 
framework of a live Galactic potential, we expect to see 
the more survival of star clusters as compared to simu¬ 
lated star clusters, which evolve in a galaxy with a con¬ 
stant mass components. Clearly, we do not claim that 
the exercise above represents a realistic effect of time- 
dependence of Galactic potential on the evolution of star 
clusters. Investigating the fate of a star cluster within a 
galaxy which grows with time and com paring wit h the dif¬ 
ferent TV-body methods (e.g. amuse ; Rieder et al.ll2013l and 
nbody6tt: iRenaud, Gieles. fe Boilyl l20114 is our upcoming 
project (Zonoozi et al., in preparation). 
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